Auxiliary particle theory of threshold singularities in photoemission and X-ray 
absorption spectra: test of a conserving T-matrix approximation 
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We calculate the exponents of the threshold singularities in the photoemission spectrum of a deep 
core hole and its X-ray absorption spectrum in the framework of a systematic many-body theory 
of slave bosons and pseudofermions (for the empty and occupied core level). In this representation, 
photoemission and X-ray absorption can be understood on the same footing; no distinction between 
orthogonality catastrophe and excitonic effects is necessary. We apply the conserving slave particle 
T-matrix approximation (CTMA), recently developed to describe both Fermi and non- Fermi liquid 
behavior systems with strong local correlations, to the X-ray problem as a test case. The numerical 
results for both photoemission and X-ray absorption are found to be in agreement with the exact 
infrared powerlaw behavior in the weak as well as in the strong coupling regions. We point out a 
close relation of the CTMA with the parquet equation approach of Nozieres et al. 



PACS numbers: 71.27.+a, 71.10.-w, 75.20.Hr 
I. INTRODUCTION 

The core level spectral function Ad(e) of a localized 
core orbital immersed in a conduction electron sea, as 
observed in the photoemission of electrons after X-ray 
absorption has long been known to show nonanalytic 
threshold behavior characterized by fractional power laws 
Ad(e) oc e~ ad in the frequency distance to the threshold 
e = uj — Eq . As shown by AndersonEI, this can be under- 
stood by considering that the sudden creation of a deep 
hole in the electronic core of an ion in a metal (or the 
filling of an empty core state) disturbs the Fermi sea of 
the conduction electrons so strongly that the subsequent 
relaxation into the new ground state follows a fractional 
power law in time rather than the usual exponential de- 
pendence. This is due to the fact that the ground states 
of the initial state and the final state are orthogonal in 
the limit of an infinite system ( "orthogonality catastro- 
phe"). At finite, but small e the relaxation process in- 
volves excitation of a large number of particle-hole pairs 
out of the Fermi sea of conduction electrons. A similar 
situation arises at the X-ray absorption threshold. There 
it has been argued that in addition to the above an. ex- 
citonic effect appears, as first discussed by MahanrJ. A 
theoretical description requires the use of infinite order 
perturbation theory. 

The problem is in some sense the simplest situation 
in which strong electron correlations are generated by 
a sudden change of electron occupations of a level cou- 
pled to a Fermi sea. The same generic problem is at the 
heart of the Kondo problem, or generally speaking, of 
quantum impurity problems, which can be understood 
as a succession of X-ray edge problems generated by suc- 
cessive flips of the impurity spin or pseudospin. In an 
even more general context, such problems arise in lattice 
models of correlated electrons, when the hopping of an 
electron from one site to the next changes the occupation 
of these sites, causing a corresponding rearrangement of 



the whole Fermi system. Given the existing evidence 
that high temperature superconductors, heavy fermion 
compounds and other metallic systems are governed by 
strong electron correlation effects, which are at present 
only poorly understood, there is an urgent need for gen- 
erally applicable theoretical methods capable of dealing 
with these complex situations. 

A powerful method of many-body physics, which di- 
rectly addresses the consequences of a change in oc- 
cupation number of a local level is the pseudoparticle 
representational! Within this framework one introduces 
pscudoparticles for each of the states of occupation of 
a given energy level, i.e. fermions for the singly occu- 
pied level and bosons for the empty level. It is well 
known that a representation of this type for the infinite 
U Anderson model of a magnetic impurity in a metal 
can give surprisingly good results already in second- 
order self-consistent perturbation theory ["non-crossing- 
approximation" (NCA-)] in the hybridization of local level 
and conduction bandu. However, at low temperatures 
and low energies the NCA fails to control the infrared 
singular behavior of the pseudoparticle spectral functions 
at threshold. Application of the NCA to the problem of 
the core hole spectral function gives a threshold expo- 
nent ad independent of the occupation of the core state, 
in contradiction with the exactly known result. 

We have recently developed an approximation Scheie 
which appears to overcome the difficulties of NCAB 
It is based on the idea of including singular behavior 
emerging in any of the two-particle channels. There 
are two relevant channels, the pseudofermion-conduction 
electron and the slave boson-conduction electron chan- 
nel. In both channels the ladder diagrams are summed, 
the resulting T-matrices are self-consistently included in 
the self-energies, as is required within a conserving ap- 
proximation scheme. The main results of this conserving 
T-matrix approximation (CTMA) are: (i) the (exactly 
known) infrared threshold exponents of the pseudopar- 
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ticle spectral functions are recoveredfl, (n) the thermo- 
dynamic quantities spin susceptibility and specific heat 
show, local Fermi liquid behavior in the single channel 
casetl and (Hi) in the -multi channel case, non-Fermi liq- 
uid behavior is foundB, in quantitative agreement with 
exact results available in certain limiting cases. 

One of the most stringent tests of a many-body method 
is the calculation of the core hole spectral function. In 
this paper we report the results of an application of the 
CTMA to this problem. 

The organization of the paper is as follows. In section 
U, we summarize the most imnoKtant results of the ex- 
act solution of the X-ray modelBEi, notably those for the 
threshold exponents for the photoemission and the X- 
ray absorption. Then, in section [II, we recall the pseu- 



doparticle representation of a spinless Anderson impurity 
Hamiltonianliil and point out its equivalence to the X-ray 
model in the infrared limit. The conserving pscudopart- 
cle approximation up to infinite order in the hybridiza- 
tion V is discussed in section [[V] and compared with the 
parquet equation approach of Nozieres et alLj in section 
The numerical results are discussed in section VI 



In 



appendix [X] we give explicitly the self-consistent equa- 
tions which determine the auxiliary particle self-energies 
within the CTMA. 



II. PHYSICAL MODEL 

The absorption of an X-ray photon by a deep level 
core electron and the subsequent emission of the electron 
leaves a core hole, which is seen by the conduction elec- 
trons as a suddenly created screened Coulomb potential. 
The simple&tjHipdel Hamiltonian describing this situation 
is given byElaJl3-y 

kcr a 

(1) 

where c ko . (cj^) are the conduction electron field opera- 
tors for momentum and spin eigenstates |k<r), with en- 
ergy £k and chemical potential /j,. The energy of the deep 
level is Ed, and V d is the screened Coulomb interaction 
between the conduction electrons at the site of the hole 
(cq ct , c 0cr ) and the hole (with operators eft, d; the spin 
state of the hole is irrelevant here). We assume that 
the hole is localized and does not have internal struc- 
ture, i.e. we neglect the finite life time of the hole due to 
Auger effect as well as a possible recoil of the hole. The 
Coulomb interaction between the conduction electrons is 
absorbed into a quasiparticle renormalization. 

Photoemission. — The spectral function of the hole, 
Ad(e), which can be measured in photoemission experi- 
ments, is obtained from the one-particle core hole Green's 
function Gd{t) = — i(T[d(t)dJ (0)]) , subjected to the ini- 
tial condition that the core hole occupation number 



d*d = for times t < (before the photoemission pro- 
cess), by taking the imaginary part of its Fourier trans- 
form, Ad(oS) — (l/7r)ImGd(w — iO). The initial condi- 
tion is equivalent to the trace (• • •) in the definiton of 
Gd(t) being taken only over states with hole occupation 
equal to zero. It is this restriction which implies the non- 
trivial dynamics of the X-ray problem. Ad(co) is propor- 
tional to the spectral weight of processes, where a pho- 
ton is absorbed by the metal, subsequently emitting the 
deep level core electron. The energy w required for this 
process is bounded from below by the threshold energy 
Eq = Ep — Score ^ AE, where E COTC and Ep are the core 
level energy and the Fermi energy, respectively, and AE 
is a renormalization due to core hole-conduction electron 
interactions. In the following we will choose the zero of 
energy such that Eq = (i.e. e = u> — Eq). The spectral 
function A d (e) then shows singular threshold behavior 
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In a landmark paper Nozieres and De Dominicisu showed 
that the exponent ad depends only on the scattering 
phase shift r\ of the conduction electrons off the core hole 
and calculated it as (s-wave-scattering) 



a d = l- (~) = l-n d , 



(3) 



where Friedel's sum rule rj = irn d has been used to ex- 
press r] in terms of the occupation number of the core 
level, rid. 

X-ray absorption. — The X-ray absorption cross sec- 
tion is given by the two particle Green's function Gi it) = 
-ie(t)([d^{t)c o<J (t),cl a {0)d(0)]} as da/de cx ImG 2 (e - 
iO). The absorption cross section is finite for e > and 
again shows singular threshold behavior 



da _ C a 
~de ~ 



(6^0+) 



(4) 



The exponent a a has been calculated by Nozieres and De 
DominicisEI with the result 



2 V [V\ 2 2 
a a = =2n d -n d 

7T \7T/ 



(5) 



III. PSEUDOPARTICLE REPRESENTATION OF 
THE X-RAY MODEL 

As will be seen below, it is useful to formulate the core 
hole problem in terms of pseudoparticles in order to im- 
pose the initial condition. We define fermion operators 
/ + (/) and boson operators b + (b) creating (annihilat- 
ing) the occupied or empty core level. The transition 
amplitude V of an electron from the core level into the 
conduction band describes the hybridization of these two 
systems. The Hamiltonian of this system takes the form 
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of an Anderson impurity Hamiltonian for spinless parti- 
cles (spin degeneracy N = 1): 



E d Pf + V ( fhc + h. c.) + XQ 



(6) 



where Co = 5Z k Ck annihilates a conduction electron at 
the impurity site. The constraint Q = pf + tfb = 1 
ensuring that the core level is either empty or occupied 
is implemented by adding the last term in (p) , where A is 
associated with the operator constraint Q = 1 and may 
be interpreted as the negative of a chemical potent^ 
the pseudoparticleso. As has been shown previously 
the limit A — > oo imposes the constraint exactly and is 
equivalent to taking all expectation values of pseudopar- 
ticle operators in the Hilbcrt subspace with Q = (no 
core hole present). Thus, in the present context, it im- 
plements exeactly the X-ray initial condition of sudden 
creation of the core hole. The auxiliary particle Green's 
functions are expressed in terms of their self-energies 



as G / 1 (iw„) = G° f b (iuj n ) 



T f,b\ 



S/(iw„), G b 1 (iv m ) 



[Gl(iv m )] - T, b (iv m ), where G° f (iuj n ) = l/(iu n - E d ) 
and G®(iv m ) — l/iv m are the respective non-interacting 
Green's functions and iu> n = (2n + l)ir/0, iv m = 2mir//3 
denote the fermionic and bosonic Matsubara frequencies. 

In the model (^]) one may distinguish two distinct 
regimes, where the impurity occupation number nd at 
infinitely long time after suddenly switching on the in- 
teraction is large (rid —* 1, Ed < 0) or small (rid — * 0, 
Ed > 0). Since, due to the hybridization, n d is equal and 
opposite in sign to the change of the conduction electron 
number (i.e. screening charge) induced by the presence of 
the impurity, rid = ~An c , these regimes correspond via 
the Friedel sum rule to large (77 — > tt) and small (rj — > 0) 
scattering phase shifts, respectively (see detailed discus- 
sion below), and may, therefore, be termed the strong 
and the weak coupling regions. We now show the formal 
equivalence between the X-ray model Eq. (Q) and the 
slave particle Hamiltonian Eq. (^) at low energies both 
in the weak and in the strong coupling regions. 

In the strong coupling region, an effective low-energy 
model is derived from the Anderson Hamiltonian (|^) 
by integrating out the slave boson degree of free- 
dom (or, equij^alently, by means of a Schrieffer- Wolff 
transformation^ onto the part of the Hilbert space in- 
volving only states with the core level occupied). The 
interaction term in the resulting effective action reads 



(7) 



x cl(iu n - iv m )c (iuj n )f(iuj n )f(iuj n + iv m ) 



where, in addition, the projection onto the physical 
Hilbcrt space is imposed by taking A — > 00. At low 
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FIG. 1. Diagrammatic representation of the effective 
low-energy interaction in the strong coupling regime, Eq. (Q) , 
and its contraction to a density-density interaction at low ex- 
citation energies. Solid, dashed and wiggly lines correspond 
here and in the following to conduction electron, pseudo- 
fermion and slave boson propagators, respectively. 

excitation energy relative to the core level, i.e. when the 
conduction electron energies after analytical continuation 
are |w|, \u' — v\ <C \Ed\ and the pseudofermions have ener- 
gies to' , lo + v w Ed (see Fig. |l|), the non-interacting slave 
boson Green's function in Eq. (0) is taken at v w Ed 
and thus reduces to 1/Ed- The resulting effective Hamil- 
tonian is thus given by Eq. ([!]), with electron operators 
eft , d replaced by pseudofermions /, interacting with 
the conduction electrons via the repulsive, instantaneous 
potential V d = -V 2 /E d > 0. 

In order to derive the effective low-energy Hamilto- 
nian in the weak coupling domain (nd — > 0, Ed > 0), it is 
useful to observe that the model Eq. (^|) is in the physi- 
cal Hilbert space invariant under the special particle-hole 
transformation / <-» b, c <-> c' and E d — > —Ed. Integrat- 
ing out the high energy states, i.e. the fermionic degrees 
of freedom in this case, and then performing this particle- 
hole transformation, the resulting low-energy Hamilto- 
nian is again given by Eq. ([l]), with the replacement 
d' , d — > /t, /, and the attractive interaction potential 
Vd = —V 2 /Ed < between conduction electrons and 
local pseudofermions. 

Having, thus, established the formal connection be- 
tween the original X-ray model Eq. ([!]) and the auxiliary 
particle Hamiltonian (|6|) in the weak and in the strong 
coupling regions, we now turn to showing that the pho- 
toemission and X-ray absorption spectra are given by the 
slave boson and the pseudofermion spectral functions, re- 
spectively. 

Photoemission. — The retarded Green's function 
Gf (t) = -ze(^([6(i),& f (0)]-) describes the propagation 
of the empty d-level in time. The corresponding spectral 
function after projection onto the physical sector Q = I, 
A b (u) = — lim^oo ImGf (<jj)/tt can be represented in 
terms of the exact eigenstates of the systemjKUhout the 
d-level, |0,n), and with the d-level, |l,n), aalS'Ll 

A b (u) = (8) 
^— £ m , n I (1, ™\b + 10, n) | 2 e-^."«5(e + e , n - e 1>m ) . 

At zero temperature (j3 = 1/T = 00), A&(e) is zero for 
e = lo — Eq < 0, where Eq = ei t Q — eo.o is the difference 
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of the ground state energies for the Q — 1 and Q = 
systems. Near the threshold, e > 0, Ab(e) has a power 
law singularity (infrared divergence), Ab(e) oc e~" 6 , for 
exactly the same reason as the hole spectral function 
Ad(e) considered above: the states \0,n) (free Fermi sea) 
and |l,n) (Fermi sea in presence of a potential scatter- 
ing center) are orthogonal, giving rise to the orthogo- 
nality catastropheo. The exponent ctb is therefore given 
in terms of the phase shift ijb (for s-wave scattering) as 
ai, = 1 — (r/i, . Using the Friedel sum rule and the fact 
that in the photoemission process (boson propagator) the 
impurity occupation number changes from initially to 
rid > in the final state, we obtain the characteristic 
dependence on nd, 



1 



(9) 



We may conclude that the threshold behavior of the phys- 
ical hole spectral function Ad(e) and the slave boson spec- 
tral function Ab(e) is governed by the same exponent, 
otd = otb, provided the scattering phase shift is the same. 

X-ray absorption. — In a similar way, the threshold 
behavior of the X-ray absorption cross section da/de 
may be obtained from the pseudofermion Green's func- 
tion. As shown in section [IJ, da/de is proportional to 
the imaginary part of the two particle Green's func- 
tion G 2 {t) = -ie{t)([dHt)c Oa (t) 7 cl a (0)d{0)}). The cor- 
responding quantity here is the slave boson-conduction 
electron correlation function 



G bc (t) = -iO(t)([b(t)c (t), 4(0)^(0)]) 



(10) 



which is given in terms of the pseudofermion Green's 
function G/(e) (after Fourier transformation) as 



G bc (e) 
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(G° f (e)y L G f (e)-l (G° f (e)y L . (11) 



It follows that the spectral functions are related by 
Abc(e) oc Af(e) ~ e _a/ , i.e. the X-ray absorption expo- 
nent is identical to the pseudofermion threshold exponent 
ctf. The latter is again determined by the orthogonality 
catastrophe argument, considering that the initial state 
of the system is now the conduction electron Fermi sea 
plus the filled d-level. The phase shift 77/, again given 
via the Friedel sum rule as the change of the occupation 
number from the initial to the final state, is now different, 
Vf — i n d — 1)tT) leading to the expression 



a f 



2n, 



d-n d 



(12) 



Comparison with (^|) again shows that the infrared be- 
havior of the pseudofermion spectral function is indeed 
identical to that of the two particle Green's function G2, 
as expected. 

It should be mentioned that in the intermediate cou- 
pling or "mixed valence" domain, ttN(Q)V 2 ks \Ed\ (rid ~ 
1/2), a Schricffcr- Wolff type projection is no longer valid 
because of large level occupancy fluctuations. The formal 



derivation of the X-ray model ([l]) from the pseudoparticle 
model (^) in the "mixed valence" regime involves a re- 
tarded effective interaction, in contrast to Eq. ([!]) . How- 
ever, since the Hamiltonian Eq. (^) is a faithful represen- 
tation of a non-interacting system (via the identification 
rft = f^b), where the constraint Q = y f+b'b = 1 merely 
serves to implement the X-ray initial condition of sud- 
den switching on the interaction between localized states 
and the conduction electrons (see above), the system is 
described by single-particle wave functions even in the 
valence fluctuation regime of this spinless model. The 
analysis of the pseudoparticle threshold exponents a&, 
a f in terms of the corresponding scattering phase shifts 
rib, rjf and the Friedel sum rule, as given above, then also 
applies in the valence fluctuation regime. It has been ver- 
ified explicitly by a numerical renormalization group cal- 
culation of the pseudoparticle threshold exponents that 
their rid dependence, given in Eqs. (|9|), (p^), is valid over 
the_complete range of the core level occupation number 
nil 

The preceding analysis shows explicitly that in the 
auxiliary particle representation the threshold exponents 
of both the X-ray photoemission and absorption are de- 
termined by the infrared behavior of single-particle prop- 
agators, involving the physics of the orthogonality cflias,- 
trophe for auxiliary bosons or pseudofermions only ! 18 !! 19 ! 
There is no separation into single particle effects and ex- 
citonic effects. 



IV. CONSERVING THEORY 



In the previous section we reformulated the core hole 
problem by introducing auxiliary particles and showed 
on general grounds that the threshold exponents of X- 
ray absorption and photoemission spectra can be ex- 
tracted from one particle properties, namely the auxiliary 
fermion and slave boson Green's functions respectively. 
In this section a systematic self-consistent approximation 
is formulated to calculate these functions. 

As a minimal requirement the constraint Q = 1 has to 
be fulfilled in any approximate theory. The constraint 
is closely related to the invariance of the system un- 
der a simultaneous local (in time) gauge transformation 
/(r) -> e e ( T )/(r), 6(r) -> e e ^b(r). The Lagrange mul- 
tiplier A assumes the role of a local gauge field and trans- 
forms as A — > A + idO/dr. Any approximate scheme 
respecting the gauge symmetry will preserve the charge 
Q in time. The simultaneous transformations /(r) — » 
e e M/(r), Ck(r) e e M Ck (r), M (r) -» M (r) + idQ/dr 
lead to the conservation of the total fermion number 
n f + J2k c k c k ~ const, where fi is the chemical potential 
of the conduction electrons (we choose ^ = 0) . Any the- 
ory which preserves these symmetries is called conserv- 
ing and may be generated by functional derivation from 
a generating functional $ of closed skeleton diagramsEl. 
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NCA. - We are interested in the limit of weak hy- 
bridization V . So let us first consider the lowest order ap- 
proximation. The conserving approximation scheme re- 
quires the self-energies to be determined self-consistcntly, 




(a) (b) (c) 

FIG. 2. (a) Diagrammatic representation of the NCA gen- 
erating functional, (b) and (c) display the pseudofermion and 
slave boson self-energies derived from the NCA functional by 
functional derivation. 

which amounts to an infinite rcsummation of perturba- 
tion theory even if only the lowest order skeleton di- 
agram ist kept (which n is known as the "non-crossing- 
approximation" (NCA)eI, see Fig. |). The NCA is known 
to yield good results in the absencer-af-,or sufficiently far 
away from a Fermi liquid fixed pointEirQ. Hence the NCA 
is not appropriate in the X-ray problem. The reason is 
that no parquet diagrams (see Fig. ||) are included in 
the lowest order approximation. By functional deriva- 
tion of $ one obtains for the slave particle self-energies 
£/ = 5&/5Gf, S;, = S$/8G b which are diagramatically 
given in Fig. and yield the set of coupled integral equa- 
tions 

/OO 7 
— G b (e + u)A c (-u)f(u) 
-oo 7T 

/OO J 
™G f (u + e)A c (u)f(u) (13) 
-oo 7T 

where A c (e) is the non- interacting local conduction elec- 
tron spectral density. At zero temperature T = the 
integral equations can be rewritten as ordinary differen- 
tial equations (with a constant density of states for the 
conduction electrons and for e — > 0)c3 

|_J_. W( 0)^ A ( f ) 



alMT) ~ »(°>^> ■ <"> 

The solution displays the well-known infrared singulari- 
ties A ftb (e) oc e~ a i- b (e -» 0) where a/ )6 = 1/2. These 
exponents obviously differ from the exact results dis- 
cussed before [Eqs. and fll^)]. 

Hence the NCA is not even in qualitative agreement 
with the exact Fermi liquid properties of the model; it 
shows no dependence of the exponents on the filling fac- 
tor rid of the deep level. This is due to the lack of vertex 
corrections which have to be included in infinite orders of 
perturbation theory, because it can be shown by power- 
counting arguments that there are,-no corrections to the 
NCA exponents in any finite orderEi 




FIG. 3. Diagramatic representation of the Bethe-Salpeter 
equations for the T- matrices in Eqs. (|l^) and (^6|), respec- 
tively. The analytically continued equations, which are calcu- 
lated numerically, are discussed in appendix |a} 

CTMA. - We have to include the major singulari- 
ties in each order of self-consistent perturbation theory. 
These singularities emerge in the conduction electron and 
pseudofermion T-matrix (Tf c ) as well as in the conduc- 
tion electron and slave boson T-matrix (T bc ). In order 
to preserve gauge invariance, self-consistency has to be 
imposed: the self-energies are functionals of the Green's 
functions which in turn are expressed in terms of self- 
energies, closing the set of self-consistent equations. The 
summation of the corresponding ladder diagrams can be 
performed by solving the integral equations, for the T- 
matrices for the pseudofermions (see Fig. phcl 



I 

T fc (iuj n , ioj' n ; itt m ) = V 2 G b (iu n + iJ n - ifl m ) 
V 2 

- -a- 2_j G b( iuJ n + i^n - iCl m )Gf(iuj^)G c (iQ m - iJ^)T Sc {iJ^ iw' n ; iQ. m ) , (15) 

and the slave-bosons 



Tbc(ii>m, iv' m \ ittn) = V 2 Gf(iv m + iv' m - iCl n ) 
V 2 x , 

JfZ^ G f( iu m + ™" n - iQ n )G b (iv'^)G c (-iD, n - iv'^)T bc (iv'^ w m \ iQ n ) . (16) 

if'' 
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Here u) n ,u}' n ,uj'^ are fermionic frequencies (w n = (2n + 
1)tt/P), v m , v' m , v'^ are bosonic frequencies [y m = 
2rmr/P), and the center of mass frequency £l m .n is 
bosonic in the case of Tf c and fermionic for T\, c . The 
self-energies £ f and £{, 



S/(iw n ) 



E 6 (w m ) = Zf CA (iv„ 



) + S{ c (^„) + Ejf (iw n ) 



7 

7c 



) + £f (w m )+S^ c (w m ) 



(17) 
(18) 



calculated from Tj c and T{, c , then follow from a gener- 
ating functional $ (see Fig. ||) by functional derivation. 
The explicit expressions are given in appendix [a|. 





FIG. 4. Diagrammatic representation of the CTMA gener- 
ating functional. The free energy diagram with two conduc- 
tion electron lines does not appear, since it is not a skeleton 
diagram. 



V. COMPARISON WITH RENORMALIZED 
PARQUET EQUATIONS 

The CTMA is closely related to the .parquet equation 
approach by Nozieres et al. In Ref. |M] these authors 
investigate the X-ray model (Q) by the methods of per- 
turbation theory. Even to the lowest order one must sum 
the so-called parquet diagrams, in close analogy with the 
Abrikosov theory of the Kojado effectcj. In this approxi- 
mation Mahan's predicitonH of the singularity in the X- 
ray absorption-spectrum was first confirmed. In a suc- 
ceeding papeilij the many-body approach was general- 
ized to include self-energy and vertex renormalization in 
a self-consistent fashion. This self-consistent formalism 



describes the reaction of divergent fluctuations on them- 
selves, and should, therefore, be useful in other more 
complicated problems, such as the Kondo effect. 





(a) 



i 1 ! 



* 2 T 



— -4 — — 




(l>) 



FIG. 5. (a) Vertex renormalization ancj-Kelf-energy repro- 
duced from the parquet equation approachu. These diagrams 
are obtained from the corresponding ones in (b) by contract- 
ing the boson-lines 1, 2 and 3. The CjEMA, therefore, con- 
tains the parquet contributions of Ref. fj] as a diagrammatic 
subclass. 



In Ref. 0] it is shown that the significant contribu- 
tions in logarithmic accuracy to the renormalized inter- 
action and the deep level self-energy are given by the 
diagrams reproduced in Fig. || (a). Both graphs are in- 
cluded in the CTMA (see Fig. | (b)): By collapsing the 
boson lines into points, i.e. by integrating out the high 
energy bosonic degree of freedom in the strong coupling 
region (n c i — ► 1) as done in section III, it is seen that the 
X-ray interaction kernel (Fig. || (a), left) can be extracted 
from the Tf, c -matrix, and the deep level self-energy (Fig. 

(a), right) is already included in the NCA. For weak 
coupling {rid ~> 0) analogous results are obtained by in- 
tegrating out the pseudofermionic degree of freedom and 
then interchanging bosons and fermions, compare section 
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[II. The self-consistent evaluation of these diagrams rep- 
resents the renormalized parquet analysis for the pseu- 
doparticles. The advantage of our formulation is that it 
is valid both in the weak coupling and in the strong cou- 
pling regime, with symmetrical expressions in these two 
regions. The symmetry between weak and strong cou- 
pling is also visible in the results for the threshold ex- 
ponents (Fig. 0). Since the CTMA is not restricted to 
parquet diagrams (which give the right asymptotic be- 
haviour only for V — > 0), but goes beyond the parquet 
approximation, one may expect that its validity extends 
beyond the weak and the strong coupling limits and in- 
terpolates correctly between these regimes. This will be 
seen the following section. 
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VI. NUMERICAL RESULTS 



The self-consistent solutions are obtained by first solv- 
ing the linear Bethe-Salpeter equations (|lf]) and (|T^) for 
the T-matrices by matrix inversion on a grid of 200 fre- 
quency points. First we insert NCA Green's functions 
into the T- matrix equations. From the T-matrices the 
auxiliary particle self-energies S/ and are calculated 
corresponding to Eqs. (Q) and (]To|), which give the re- 
spective Green's functions. This process is iterated until 




.4 I 1 1 1 1 1 1 I 

-6 -5 -4 -3 -2 -1 
lg e 



FIG. 6. Auxiliary particle spectral functions Af and Ab in 
the weak coupling regime in a logarithmic plot. The ener- 
gies are in units of the half band width D. The slopes of the 
dashed lines indicate the exact threshold exponents. 

convergence is reachedB. The T-matrices show nona- 
lytic behavior in the infrared limit. 

As can be seen from Fig. || the fermion and boson 
spectral functions display power law behaviour at low 
frequenciesti The power law behavior emerges in the 
infrared limit, i.e. for energies smaller than the low en- 
ergy scale (which is Ed)- For smaller frequencies there 
is always a deviation from the power law behaviour due 
to finite temperature. The exponents extracted from the 
spectral functions at low but finite temperature for var- 
ious values of the deep level filling rid in Fig. [7] are in 
good numerical agreement with the exact results in the 
regions rid & [0.0,0.3] and rid S [0.7,1.0]. Note that 
in contrast to the n^-dependent exponents within the 
CTMA the NCA spectral functions always diverge with 
n^-independent exponents a/ = ab = 1/2. For inter- 
mediate coupling, n d <E [0.3,0.7], the convergence of the 
self-consistent scheme is very slow, and we find no sta- 
ble numerical solution. It remains to be seen whether 
this is due to numerical instabilities or possibly due to 
the importance of further vertex corrections beyond the 
CTMA. 




0.2 0.4 0.6 0.8 1 

FIG. 7. Auxiliary particle threshold exponents exctracted 
from spectra as in Fig. [sjfor a number of deep level fillings rid- 
The solid lines represent the exact values derived in Eqs. (^) 
and 

A comparison of the CTMA results with the weak- 
coupling treatment, which corresponds to rid — * in our 
model, shows that for finite interaction strength renor- 
malization effects are important (see Fig. |^). The con- 
nection between rid and Ed/T is exactly given by Friedel's 
sum rule rid — 1/2 — arctan(£^/r)/-7r. Again we mention 
the rid dependence of the exponent af in contrast to the 
NCA result: To recover the Fermi liquid properties of 
the model one thus has to go far beyond the lowest order 
self-consistent approximation. 




I 1 1 1 1 1 

3 2.5 2 1.5 1 

rVr=1/[7iN(0)V d ] 

FIG. 8. Comparison psf,— the CTMA results and the 
weak-coupling calculationtixta for the threshold exponent of 
X-ray absorption spectra. 



VII. CONCLUSION 

In summary, we have calculated the exponents of 
threshold singularities in the X-ray photoemission and 
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absorption spectra, using a standard many-body tech- 
nique, where the empty and the singly occupied core level 
are represented by separate fields, auxiliary bosons and 
pseudofcrmions, respectively, coupled to the conduction 
electrons via a hybridization interaction. In this formu- 
lation, the X-ray problem is described by a spinless An- 
derson impurity model in pscudoparticle representation, 
and the initial condition of sudden creation of the impu- 
rity potential is implemented by the constraint that all 
expectation values of local fermion or boson fields must 
be calculated in the Hilbert subspace with pseudoparticle 
number Q = 0. The latter can be fulfilled exactly. It was 
further shown that the X-ray photoemission cross section 
or core level spectral function is given by the boson spec- 
tral function, while the X-ray absorption cross section 
is proportional to the total fermion hybridization ver- 
tex. Therefore, the X-ray photoemission and absorption 
threshold exponents are identical to the infrared expo- 
nents of the auxiliary boson and pseudofermion spectral 
functions, respectively. It follows that both X-ray pho- 
toemission and absorption are solely governed by the or- 
thogonality catastrophe, and there is no separation into 
single particle and excitonic effects. 

In a more general context, the generalized 
SU(iV) xSU(M) Anderson impurity models, classified by 
the spin degeneracy N of the local orbital and the number 
M of degenerate conduction electron channels, may be 
considered as standard models to describe strong correla- 
tions induced by the restriction of no double occupancy 
of sites. Depending on their symmetry, these models 
display Fermi (N = M — 1 or N > M + 1) or non=Fermi 
liquid behavior (2 < N < M) at low temperatureEj. The 
present case of the spinless Anderson impurity model in 
slave boson representation (N — 1, M = 1), Eq. (||), 
may be considered as the most stringent test case for 
the development of new methods for strongly correlated 
systems. This is because for this case earlier approxima- 
tion schemes like the non-crossing approximation (NCA) 
fail in the most pronounced way to even qualitatively 
describe the low-energy Fermi liquid behavior of this 



model, i.e. the rid dependence of the infrared threshold 
exponents, while in the non-Fermi liquid case the NCA 
gives the correct, exponents at least in the Kondo limit 
of these modelsEi 

In the present paper we have applied a recently de- 
veloped approximation scheme, the conserving T-matrix 
approximation (CTMA) to the N = 1, M = 1 Ander- 
son impurity model to calculate the X-ray photoemis- 
sion and absorption threshold exponents on a common 
footing. The CTMA includes the complete subclass of 
diagrammatic contributions which, in the limits of weak 
(rid ~ * 0) and strong (rid 1) impurity scattering poten- 
tial, reduce to the renormalized pa*ouet diagrams, which 
have been shown by Nozieres et alJl3 to describe the exact 
infrared singular behavior in the weak coupling regime of 
the X-ray problem. As a result, the CTMA recovers the 
correct X-ray photoemission and absorption exponents 
in a wide region around weak as well as strong coupling. 
In connection with earlier resultsQ on the spin 1/2 An- 
derson impurity model (N = 2, M = 1), this makes 
the CTMA the first standard many-body technique to 
correctly describe the Fermi liquid regime of the Ander- 
son impurity models in a systematic way, including the 
smooth crossover to the high temperature behavior. 

We are grateful for discussions with J. Brinkmann, T. 
A. Costi and T. Kopp. T.S. acknowledges the support of 
the DFG-Graduiertenkolleg "Kollektive Phanomene im 
Festkorper" . This work was supported in part by SFB 
195 of the Deutsche Forschungsgemeinschaft. Computer 
support was provided by the John-von-Neumann Insti- 
tute for Computing, Julich. 

APPENDIX A: CTMA EQUATIONS 

In this appendix we give explicitly the self-consistent 
equations which determine the auxiliary particle self- 
energies within the CTMA. In the Matsubara represen- 
tation the vertex functions Tf c and T bc are given by the 
following Bethe-Salpeter equations: 



Tf c (iuj n ,iuj' n ;ifl m ) = If c (iw n ,iuj' n ;iQ. m ) 
V 2 x 

+ — G b (iu n + iu n - i£l m )Gf(iu n )G c (iQ, m - iu n )Tf c (itJ n , iu n \ ift m ) (1) 



with 

I fc (iuJn, iu' n \ ifim) = — o" ^2 G b( iuj n + iu'n - i&m)G f(iu'£)G c (iSl m - iu")G b (iuj' n + iu'Jl - iSl m ) , 



and T bc 

T bc (iv m ,iv' m ;iQ, n ) = I bc (iv m ,iv' m ;iVl n ) 
V 2 

o-yi G f(iv m + iv" n - £O n )G 6 (u4)G c (iz/^ - iVL n )T bc (iv'^,w l m ;iQ, n ) (2) 

iv 1 ' 
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with 



I bc {iv m ,iv' m ; iSl n ) = — — ^2 G f{iv m + iv'm ~ ^«) G fc(* jy m) G c(^m - iSl n )G f (ii/ m + iv" n - ift n ) . 

iv 1 ' 



Note that, in addition to the different sign in Tf c , these with two bosonic or fermionic rungs, respectively. The 
vertex functions differ from the T-matriccs defined before terms with a single rung correspond to the NCA diagrams 
in that they contain only terms with two or more rungs, and are evaluated separately, 
since the inhomogenous parts If c and I bc represent terms The fermion self-energies in Fig. ^ are given by 
I 



Yj f f c (ioj n ) = — ^ G c (iSl m - iuj n )Tf c (iu; n ,iu> n ;iQ m ) (3) 
V 2 

Y, b f c (iuj n ) = -— ^2 G c (iu n - ii/' m )G b (ii/ m )G c (iuj n - iv^)G b (iv^ (4) 

and the boson self-energies by 

%b c (ivm) = 4 X! G c {iv m -iSl n )T bc {iv m ,iv m ;iSl n ) (5) 
V 2 

Z f b c (iuj n ) = -— G c(iu' n -iv' m )Gf{iw' n )G c {iu'l-iv m )Gf{iw'^ (6) 



After analytical continuation to the real frequency axis we have to solve the NCA equations (13) and the following 
CTMA equations 



Ej°(e)=/ ^f(u-e)A c (u-e)T fc (e,e;u) (7) 



S^ c (e) = -V 2 I — —f(u- e)f(u' - e)A c {e - u')G b (e)T bc (u, u'; u + u' - e)A c (e - u')G b {u') (8) 



OO J — oo 



/CO J 
™f( u -e)A c (e-u)T bc (e,e;u) (9) 
-oo T 

/oo j f oo j i 
— —f(u- e)f(u' - e)A c (e - u)G f (e)T fc (u, «'; u + u'- e)A c {u' - e)G f (u') (10) 
-oo T J-oo T 

with the fermion-conduction electron vertex function 

/oo 1 
— f(u - Sl)G b (e + u- Q)G f (u)A c (Q - u)T fc (u, e'; SI) (11) 
-oo K 

/OO J 
— f(u - Sl)G b {e +u- Sl)G f {u)A c {Sl - u)G b {e' +U — SY) 
-oo 7T 

and the boson-conduction electron vertex function 



/OO j 
— f(u - Sl)G f (e + u- Sl)G b (u)A c (u - Sl)T bc (u, e'; fi) 
- oo 7T 

/CO 7 
— f(u - Sl)G f (e + u- Sl)G b {x)A c {u - Sl)G s {e' + u-Q). 



(12) 







Note that the self-energy contributions obtained from the 
two rung T-matrix terms (J/ c and J{, c ) display no skeleton 
diagrams; they are subtracted in the end. 




Tfc 



FIG. 9. Diagrammatic representation of the NCA and 
CTMA expressions for the pseudoparticle self-energies. 
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